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ABSTRACT 

We present a study on Spectral Energy Distributions, Morphologies, and star for- 
mation for an IRAC-selected extremely red object sample in the GOODS Chandra 
Deep Field-South. This work was enabled by new HSTfWFC3 near-IR imaging from 
the CANDELS survey as well as the deepest available X-ray data from Chandra 4 
Ms observations. This sample consists of 133 objects with the 3.6/im limiting magni- 
tude of [3.6] = 21.5, and is approximately complete for galaxies with M* >1O*'M0 at 
1.5<z<2.5. We classify this sample into two types, quiescent and star-forming galax- 
ies, in the observed infrared color-color ([3.6]-[24] vs K-[3.6]) diagram. The further 
morphological study of this sample show a consistent result with the observed color 
classification. The classified quiescent galaxies are bulge dominated and star-forming 
galaxies in the sample have disk or irregular morphologies. Our observed infrared 
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color classification is also consistent with the rest-frame color (U-V vs V-J) classifi- 
cation. We also found that quiescent and star-forming galaxies are well separated in 
the nonparametric morphology parameter (Gini vs M20) diagram measuring their con- 
centration and dumpiness: quiescent galaxies have Gini coefficient higher than 0.58 
and star forming galaxies have Gini coefficient lower that 0.58. We argue that the star 
formation quenching process must lead to or be accompanied by the increasing galaxy 
concentration. One prominent morphological feature of this sample is that disks are 
commonly seen in this massive galaxy sample at 1.5<z<2.5: 30% of quiescent galax- 
ies and 70% of star forming galaxies with M* > lO^^M© have disks in their rest-frame 
optical morphologies. The prevalence of these extended, relatively undisturbed disks 
challenges the merging scenario as the main mode of massive galaxy formation. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high- 
redshift — galaxies: structure 



Introduction 



Understanding when and how the most massive galaxies formed remains one of the most 
outstanding problems in galaxy forma tion. In the presen t-day universe, most massive galaxies 
(M > IO^'Mq) are early-type galaxies ([Baldry et al.ll2004|) . which are primarily composed of old 
stellar populations with no or little star formation. Spectroscopic studies of large local samples 
reveal that the bul k of stellar mass in massive early- t ype galaxies was fo rmed at z > 1.5 within 
a short timescale ([Thomas et al.l 120051: lRenzinill2006l : I Jimenez et al.ll2007L and refs. therein). At 
intermediate redshifts, several surveys reveal that the precursors of present-day massive early-type 
galaxies have already turne d passive by z 0.8, but may continue to assemble the i r stellar mass 



later t hrough dry mergers (IBell et al 



20041 Yamada et al. 20051: iBundv et al.ll2006l: ICimatti et al 



20061 : iBorch et allbood: Ipaber et al.lboovl: [Brown et al.lboovl) . The lack of significant star for 



mation in these galaxies since z = 0.8 implies that their stellar mass formed at higher redshifts. 
Recently, both observational and theoretical results suggest that the cr i tical epoch for rnassive 



galaxy evolution and mass assembly is at 1 < z < 3 (IFontana et al.ll2004l: 



Nagami ne et al. 



2005 



2007; Ilbert et al. 



201o|y ^ Thus, a complete census of a well-defined massive galaxy sample at this 



Kong et al.ll2006l : |Papovich et al 



2006 



Glazebrook et al. 



2004 



Abraham et al.i,2007: Amouts et al. 



redshift range is crucial to understanding the formation and evolution of these objects. 

Significant progress has been made in detecting and identifying massive galaxies at 1 < z < 3 
using large-format near-infrared (NIR) array cameras on 4-10 m class telescopes. Recent wide- 
field NIR surveys reveal several "new" galaxy populations, which includes Extremely Red Objects 
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(EROs; R-K> S. lElston et al]ll988l) . IRAC-selected EROs ( zs^n - r3.61 >3.2fl |^n_eLalJl2QQ3), 
BzKs (sBzK:(z-K)-(B -z) > -0.2, pBzK : z-K > 2.5: lDaddi et aDl2004l) . and Distant Red Galax- 
ies (DRGs, J-K > 2.3. iFranx et al.ll2003|) . These galaxy populations are mainly selected on the 



basis of their red rest-frame optical colors, and have typical stellar masses in excess of IO^^Mq 
down to current NIR limits. Moreover, various studies show that they domin a te the high-mass 
regime of the z ~ 2 stellar mass function (SMF). For instance, IConselice et alj (l2008h suggested 
that EROs contribute 75% of all > W^Mq galaxies at z ~ 1.5 down to Kyega = 19.7, while 
van Dokkum et al.l (|2006h found a similar fraction for DRGs in the range of 2 < z < 3. This makes 
red color selection an efficient and effective way of studying massive high-redshift galaxies. 

Utilizing these novel techniques for identifying high-redshift massive galaxies, a great deal of 
work has been carried out to explore their nature. These galaxies turn to be very different from their 
local counterparts. First, star formation activities were more intense at higher redshifts. At z ~ 2, 
more than half of the massive galaxies are rapidly star-forming galaxies wit h > lOOMfr^yr"^ while 
the remaining are quiescent galaxies with no or little ong oing star- formation ( Papovich et al.ll2006 : 
Daddi et al.ll2007l : lcassata et al.ll2008l : iKriek et aLlboOsb . In contrast, at z ~ most of the galaxies 
with similar mass are quiescent galaxies. Second, while the cosmic mass density in massive star- 
forming galaxies is flat or decreases from z ~ 2 to the present day, the mass density in massive 
quiescent galaxies increases by a factor of ~ 10 over the same time period. Such dramatic evolution 
can not be fully explained by modest mass growth of individual quiescent galaxies. Transformation 
of the star-forming galaxies into quiescent galaxies is indispensable to match this d ensity evolution 



of qu iescent galaxies, as well as that of star- forming galaxies themselves (see, e.g. jBrammer et al. 



20111). Understanding this transformation is now one of the central questions in studying galaxy 
evolution. 

To study the physical mechanisms responsible for this transformation, it is essential to clas- 
sify galaxies into different types. Though a clear bimodal color distribution between quiescent 
and star- forming gal axies persists to at least z ^ 1, the birnodality becomes less distinct at the 
high-mass end (e .g., Istrateva etaPboOll: leianton et aDl2003l : baldrv et al.ll2004l : beU et all 12004 : 
Faber et al.ll2007|) . This may be largely attribut ed to an increase in dust content in the massiv e star- 
forming galaxies which reddens their colors ([Baldry et al.ll2004l l2006l : Driver et al.ll2007|) . The 
dust reddening turns to be more significant at higher redshifts due to increasing importance of 
obscured star-formation. At z ~ 2 most massive star-forming galaxies have very red optical-to- 
infrared colors. Indee d, these star-forming galaxies have such red colors that they are prematurely 
on the red sequence (IB rammer et al.ll2009|) . A number of other methods w ere proposed to sepa- 
rate the two populations, using for example, observed colors I- J vs. J-K (IPozzetti & Mannucci 



'z850 is the F850LP band magnitude. 
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colors (|Papovich et alj|2006l : iFontana et al. 



2000) . rest-frame colors U -V vs. V-J (Wuvts et all 120071: IWilliams et alJ l2009h. mid-infrared 



20091), or directly us ing output spectral types and spe 



cific star-formation rate s (SSFRcb from th e SED fitting process (|Amouts et al.ll2007l: iDamen et al 



20091 : lllbert et al.ll2010l: IGuo et alj|2011|) . Each of these techniques has its advantages and dis- 
advantages, and they can give different classification results for a significant number of galaxies. 
Resolving these discrepancies is important for us to have a clean separation of quiescent and star- 
forming populations. Moreover, the fract ion of AGNs among massive galaxies at z ~ 2 is high, up 
to ~ 30% (see, e.g.. iMessias et al.ll2010l) . AGN activity can also contribute to the infrared emis- 
sion, which is used to identify star-forming galaxies in most of the classification methods. Thus a 
comprehensive identification of AGN population is necessary for a better understanding of galaxy 
classification. 

In addition to the star-formation status, morphologies provide important constrai ns on the 



galax y assembly history and are crucial to distinguish among models of galaxy formation (|Roberts & Haynes 
994b. Studies of morphologies for high-redshift galaxies have greatly advanced in the last decade 



(IPapovich et al. 



2005 



Forster Schreiber et al. 



.otz et al. 



2011 



2006; 



Toft et al 



20071 : iBuitrago et al.ll2008l: ICameron etal.ll2011 



Kartaltepe et al.ll201 Several recent work present evidence that 
galaxy morphologies are correlated with star-forming status for galaxies at z ~ 2: quiescent galax- 



ies tend to be spheroidal a nd compact while star-forming ga 



lar/mergers with larger si zes doa ddi et aL 



2005 



2008; 


Krieketal.ll2009a: 


2011; 


Wuvts et al. 


2011; 



Damianov et al 



20091 : 



Toft et al.ll2007 



Szomoru et al 



axies are mostly disks or irregu 



van Dokkum et al. 2008 



2010l:lCassata etal.ll2011 



Buitrago et al. 



Weinzirl et al. 



Bell et al.ll201l[) . However, this correlation at z ~ 2 shows large scatter 



and is far from clear. A significant number of systems are clas sified as having the sa me star- 



formation status yet with different mo rphological types (see, e.g ., IConselice et al.ll201 la|) . Based 



on the new HST/WFC3 NIK imaging, Ivan der Wei et all (120 



Cameron et al. 



11) sh ows that ~ 65% of quiescent 



20111) . though their sample is small 



galaxies are actually disk-dominated (also see, 
. Moreover, many of these work characterized galaxy structures with the Sersic index obtained by 
fitting a Sersic model to the surface brightness data, i.e., a model-dependent approach. Yet there 
are con cerns whether the pr ofile of a high-redshift galaxy can be well fitted with a single sersic 
profile (|Szomoru et al.ll2010|) . Thus other methods, including both visual classification and model- 
independent measurement of galaxy morphology, are needed to probe the relationship between 
galaxy morphology and their physical properties. Nevertheless, this requires high-resolution and 
high signal-to-nosie NIR (rest-frame optical) imaging for a statistical sample of z ~ 2 galaxies, 
which now becomes available thanks to the Cosmic Assembly Near-infrared Deep Extragalactic 
Legacy Survey (CAND ELS; PI: Sandra M. Faber and Henry C. Ferguson, cf, iGrogin et al.ll2011 ; 
Koekemoer et al.ll201 1 ). 



2 SFR per stellar mass, SSFR = SFR/M, 
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In this paper, we present a study on star-formation activity, morphologies, and rest- frame col- 
ors for a sample of massive galaxies at z ~ 2 selected in the GOODS-South. This work was enabled 
by the new HST fWFC3 near-IR imaging from the CANDELS as well as the deepest available X- 
ray data from Chandra 4 Ms observations. These new data, supplemented with the available rich 
data set in this field, allows us to have a better separation between quiescent galaxies, sta r-forming 



galax ies, and AGNs. We use a sample of bright IRAC-sele cted Extremely Red Obj ects (lYan et al. 



20041) . which effectively reject low-redshift contaminates. iConselice et al.l (I2011bl) argue that the 
lERO selection is the most efficient among other methods (EROs, BzKs, and DRGs) for identifying 
massive galaxy samples at 1.5 < z < 3. 

The structure of this paper is as follows. The sample selection is described in section 2. We 
discuss the classification between quiescent and star- forming galaxies based on their star- formation 
status in section 3. Both visual classification and nonparametric approaches are performed on 
their H ST /WFC3 F160W morphologies in section 4. After separating quiescent and star-forming 
galaxies, we derive their redshift and stellar mass in section 5. We examine the joint distribution 
of star formation activity, morphology and rest-frame colors in section 6. A summary of our main 
results is given in section 7. Throughout the paper, we assume il\ = 0.7, Qm = 0.3, and Hq = 70 
km s"^ Mpc"^ All magnitudes are in the AB system unless specified otherwise, and the notation 
"[3.6]" means the AB magnitude at wavelength 3.6 fxm. 



2. Data Sets and lERO Selection 



The GOODS Chandra Deep Field-South CGOODS-S. iGiavalisco et al.ll2004h is among the 
best studied deep fields, with a rich multi- wavelength data set. The ultradeep IRAC imaging of 
GOODS-S covers the whole 10' x 15' region (Dickinson et al., in preparation) and has an exposure 
time of 23 hrs per pointing, yielding a 5a limiting magnitude of [3.6] ~ 25.6. Zheng et al. (2012, in 
preparation) reprocessed all GTO and legacy IRAC survey data in five CANDELS fields including 
the GOODS-S to produce a homogeneous catalogue. We adopted their 3.6 fim selected sample as 
our primary catalogue. 



Giavalisco et al.l(|2004|) performed very deep ACS imaging including F435W, F 606W, F775W, 



and F8 50LP bands in the GOODS-S, with a 5a limiting magnitude of zsso ~ 27.3. iRetzlaff et al. 



(|2010l) also carried out deep VLT/ISAAC NIR imaging in this field, with limiting magnitudes J = 
25.0, H = 24.5, and = 24.4 (5a). The deep NIR photometry i s essential for der iving accurate 
photometric redshifts and stellar masses for red galaxies at z ~ 2 dWuyts et al.ll2007h . 



MIPS 24 u rn dRieke et a 



range < z < 3 (|Huang et al 



20051 



■1120041) is very sensitive to dust emission from galaxies in the redshif t 



20091 : lYan et al.ll2005l : IWebb et al.ll2006l : IPapovich et al.ll2007h . 
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and is thus extremely useful to identify the dusty star-forming population in red galaxies. The 24 
/im imaging in the GOODS-S is the deepest ever obtained with an exposure time of 10 hour 
per pointing. The 3cr point-source sensitivity limit is ~10//Jy, permitting detection of luminous 
infrared galaxies (LIRGs) with Ljr 2x IO^'Lq out to z ~ 2.5. On the other hand, the 24 /im 
emission from a galaxy can also be powered by an AGN, which can be identified by deep X-ray 
observations. The recent 4 Ms Chandra imaging in the GOODS-S, the deepest Chandra survey 
ever obtained, reaches on-axis flux limits of ~ 3.2 x 10"'^, 9.1 x 10"'*^, and 5.5 x 10"' ^ erg cm~ ^ 



s"^ fo r the fuU (0.5-8 keV), soft (0.5-2 keV), and hard bands (2-8 keV), respectively (|Xue et al. 
2011|) . These deep X-ray data permit detection of even highly obscured AGNs. 



A major breakthrough in deep imaging of GO ODS-S comes from t\y o recent HSTfWFC3 
NIR surveys: the WFC3 Early Release Science (ERS. IWindhorst et al.ll201 11) program and the first 
phase of the CANDELS survey. The ERS mosaic covers the northern 10' x 4' region of GOODS-S 
while the CANDELS dataset covers the southern 10' x 12' region. Both surveys reach 5a limiting 
magnitude of > 27.0 in F160W (//-band), providing high-quality structural information (concen- 
trations, Gini/Min) for galaxies as faint as H = 24.6, which is equivalent to 2 x IO^^Mq at z ~ 2 
dGrogin et al.lboul) . 



The sample is selected within the central 138 arcmin^ area in the GOODS-S where deep ACS, 
WFC3, ISAAC, and IRAC data are all available. We set criteria 



[z85o] - [3.6] > 3.25 and [3.6] < 21.5 
to select a bright lERO sample. Figure [T] shows that the color crite rion 



selects quiescent galaxies and dusty star-forming galaxies at z > 1 (|Yan 



(1) 

Z8.5n] -[3.6] > 3.25 mainly 



2008). Their red colors are 



due to either the Balmer/4000 A break shifting beyond F850LP or their steep slope of the rest-frame 
UV spectra. Young and blue galaxies can not reach this threshold even at higher redshifts. Thus 
simple z85()-[3.6] cut are very effective in selecting galaxy samples at high-redshift. Moreover, the 
bright magnitude cut of [3.6] < 21.5 ensures that all galaxies in this sample have high-quality SEDs 
as well as reliable morphological measurements. Our final lERO sample includes 133 galaxies. 
The median z»5o band magnitude in our sample is zsso ~ 24.6, and 80% of them have at least a 3a 



detection in all the ACS and IRAC bands. For comparison, the lERO sample in lYan et al.l (120041) 
is much fainter with a m edian zsso ~ 27. We then matched these sources to the FIREWORKS 
catalog ( Wuyts et al.ll2008 ) and obtained multi- wavelength photometry from ground-based U band 
to Spitzer/MlPS 24 fim. 29% of the sample, or 38 sources, have an X-ray counterpart in the 4 Ms 
Chandra source catalog within a search radius of 1.5". 

Figure [T] also shows the distribution of their [3.6] - [4.5] colors for this sample, which can 
be used to constrain their redshifts. Most galaxies in the sample have [3.6] - [4.5] > 0, indicating 
that they are mainly at z > 1.5 when the 1.6 yum stellar emission bump begins to move into the 
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IRAC 4.5 band (ISawickill2002l : iHuang et alJl2004l : |Papovichll2008l : iHuang et al.ll2009h . Their 
photometric redshifts are presented in section |5] 



Separating lEROs into Quiescent and Star-forming Populations based on [3.6] - [24] 

color 



We have shown that the lEROs tend to be either quiescent galaxies (quiescents) or dusty 
star-forming galaxies (dSFGs) at z > 1.5. However, it is difficult to separate them solely based 
on observed UV-to-NIR SEDs. Instead, the mid-infrared, e.g., the MIPS 24 fim, traces star- 
forming activity across a wide range of redshifts, and can be used to separate these two populations 
dPapovich et allboodliontana et al.ll2009l : lFang et allboool) . In particular, in redshift 1 .5 < z < 2.5 
the 7.7 fim PAH emission shifts in to the observed 24 /xrn , allowing a universal conversion of 2 4 /xm 
flux into star-formation rate (SFR jPapovich et al.ll2007[ iRieke et ahlboogl : Elbaz et al.ll201ll) . On 
the other hand, the IRAC 3.6 /im i n this redshift range probes rest- frame N IR, which can be treated 
as a good proxy for stellar mass (iBell & de Jongll2001t ICole et al .1120011) . Thus we argue that the 
[3.6] - [24] color provides a good indicator of specific star-formation rates (SSFR) for galaxies in 
this sample. 

We plot K - [3.6] vs. [3.6] - [24] color-color diagram for X-ray undetected and X-ray detected 
lEROs in Figure [2^ and Figure |2j5, respectively. The lEROs show a clearly bimodal distribution 
in the [3.6] - [24] color. The 3a limiting flux density at 24 um in the GOODS-S is ~ 11/iJy, 
equivalent to a SFR of ~ 1 Mfr, yr"' at z ~ 2 dRieke et allboool) . lEROs have a typical stellar 
mass of ~ IO^^Mq (see.e.g. lYan et al.l 120041) . thus this 24 /im detection limit is equivalent to a 
SSFR ~ 10"^° yr"i. All the 24 /xm-undetected galaxies are in the [3.6] — [24] < 0.3 region in 
Figure 121 and we propose an empirical color criterion, [3.6] - [24] = 0. 3, to classify this sample 
into quiescents and dSFGs. This criterion is consistent with that used in lDamen et al. (l2009.i^ We 
also verified that by changing the [3.6] - [24] color criterion within a scatter of 0.2, the numbers 
in the quiescents and dSFGs population changed at most ~10%, and our main conclusions drawn 
in this paper remain valid. 

In several cases, 24 fim emission may be very weak or absent even for dSFGs. At z ~1.4, 
the silicate absorption at 9.1 fim shifts in the MIPS 24 fim band. T hus SFGs with a stro ng silicate 
absorption feature may have no 24 fim detection at this redshift (iMagdis et al.ll201lh . We thus 
propose another K - [3.6] color to At z ~ 2, K - [3.6] probes rest-frame R-J, which is sensitive to 



Damen et alj j2009l) adopted SSFR < 1 /(3 x f//) yr ' to select quiescents, where tn is the age of the universe at 



a given redshift. This criterion is used to characterize the low-SSFR peak of massive galaxies, which increases with 
redshift. At z = 2 this yields SSFR < lO"'" yr"' . 
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dust extinction. DSFGs have redder K - [3.6] colors compared to quiescents at similar redshifts, 
as shown by the color tracks of a quiescent and a dSFG SED template in Figure |2k. Most objects 
in our sample with [3.6] - [24] < 0.3 have blue K - [3.6] , but a few of them have very red K - 
[3.6] (see galaxies to left of dashed line in Figure |2t a)). Photometric redshifts for these sources 
confirm that they are at z ~ 1.4 or at very high redshift, z ~ 3, where 24/im is not sensitive to star 
formation. We argue that these sources are really dSFGs. We therefore propose the observable 
color criteria for identification of quiescents in the sample as: 



[3. 6]- [24] < 0.3 and i^- [3.6] < 1.2. 



(2) 



The remaining objects are denoted as dSFGs. We note that though the adopted ^- [3.6] color 
criterion is empirical, the number of these 24 /im-faint dSFGs is very small and have redshifts 
beyond our main interests. 

Based on Eq[2l a total of 90 galaxies, 68% of the whole sample, are classified as dSFGs; 
43 galaxies, 32% of the whole sample, are quiescents. The median [3.6] - [24] colors for the 
two populations are ~ 2.2 and -0.5, respectively, indicating that their SSFRs are significantly 
different. In fact, those galaxies with [3.6] - [24] > 0.3 have a median 24 jim flux density of ~ 
125 /iJy, implying that a significant fraction of them are ultra-luminous infrared galaxies (ULIRGs) 
at z ~ 2. Roughly 30% of this lERO sample, 38 sources, ar e detected at X-r ay, all of which are 
classified as AGNs according to their X-ray luminosities in IXue et al.l (|201lh . Among these X- 
ray lEROs, 79% are dSFGs and the remaining 21% are quiescents. We d o note that qui escent 
galaxies with AGNs, which can also have strong 24 /im emission (see, e.g.. IGu et al.ll2007|) . may 
be misidentified as dSFGs. As shown in Figure |2j5, several AGNs even show very blue K - [3.6] 
contrary to their red [3.6] - [24]. Further classification of these X-ray lEROs will be performed 
using their morphologies and rest-frame colors. 



4. Morphologies of lEROs from //5rAVFC3 NIR Imaging 



Galaxy morphologies provide a dditional information on the mass assembly history of galaxies 
(|Bell et al.ll2005l : IZheng et al.ll2004h . We performed both visual classifications and nonparametric 
morphological measurements to see whether quiescents and dSFGs in our sample can be separated 
based on their morphologies. All morphological analysis in this work are based on the CANDELS 
6-epoch v2.0 drizzled mosaics for the WFC3 F160W band, with a pixel scale of 0.06", 
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4.1. Visual Classification 

The lEROs in our sample show very diversified morphologies, covering a wide range of types 
from extended and disturbed low surface brightness features to bulge-dominated morphology. We 
divide morphological types for this sample into three broad categories: 

Spheroid: Single, round, and centrally concentrated source with no evidence of extended low 
surface brightness features. 

Disk: Undisturbed source with extended low surface brightness features. 

Irregular/Merger: Single highly irregular galaxy with evidence of non-axisymmetric, ex- 
tended low surface-brightness features or two or more distinct galaxies showing distortions and 
interaction features such as tidal arms. 

Figure |3] shows typical examples of the three types in the WFC3 F160W images. Visual in- 
spections for this sample were performed by three of the co-authors independently. We combined 
the classifications from each inspector and reviewed the sources together to resolve the disagree- 
ments, which involved ~ 15% of the sample. The final visual morphological types are listed in 
Table [TJ About 40% of this sample are classified as disk galaxies, 30% are spheroids and the re- 
maining 30% are irregular/mergers. The X-ray sources have the similar fractions of disk, spheroid 
and irregular/merger types: 45%, 20% and 35% respectively. Mos t of the X-ray sources show a 



bulge component even they are identified as disk galaxies (see also, iKocevski et alj|201 ID . 

We present F160W images for quiescents and dSFGs in Figures |4]and[5l respectively. Roughly 
70% of the quiescent galaxies are spheroids, and 30% of them show a significant disk compo- 
nent. The dSFGs, in contrast, are either regular disks or irregular/mergers, with only a few being 
spheroids. It should be noted that he quiescent disks are significantly different from the star- 
forming disks: the former also have a prominent bulge while the latter have small or no bulges. 
Thus we can separate the two populations based on whether they have a significant bulge compo- 
nent: Both spheroids and those disks with a prominent bulge are quiescents while the remaining 
galaxies are dSFGs. 



4.2. Nonparametric Morphological Parameters 

Nonparametric morphological parameters, such as Gini and M20, have now become popular 
in the automated identification of galaxy morphology. These parameters do not rely on certain 
model parameter fits, such as t he Sersic index, a nd therefore can be applied to irregulars as well as 
standard Hubble-type galaxies (|Lotz et al. Il2004h . The Gini coefficient is a concentration parameter. 



measuring the relative distribution of the galaxy pixel flux values. It is high for galaxies with much 
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of their light concentrated in a small number of pixels, regardless whether those pixels are in the 
projected center. Therefore, the Gini method has the advantage that it do not require the galaxy to 
be circularly symmetric, and is particularly useful for galaxies at high redshifts, which can have 
very distorted morphologies. M20 describes the second-order moment of the brightest 20% of the 
galaxy's flux, which is very sensitive to merger signatures such as multiple nuclei, tidal tails, and 
off-center star clusters. Gini and M20 are shown to be closely related with visual morphologies for 
galaxies both at low and intermediate redshifts: early-type galaxies have higher Gini and lower 
while late-type galaxies have lowe r Gini and higher M 20 (ILotz et al.l 120061 : [Abraham et al.l 12007 : 
Capak et allboovl : Lotz et allboosl : iKong et alJ 120091) . At higher redshift, however, it is difficuk 
to measure Gini and M20 due to a lack of high-quality NIR (rest-frame optical) images. And their 
relationships with visual morphologies, as well as star-forming status remain unclear. The new 
WFC3 data from CANDELS provides NIR imaging with sufficient angular resolution and depth, 
allowing us to reliably measure Gini and M20 for a representative sample of massive galaxies at 
z ~ 2. 



We calculated Gini and M20 for galaxies in our sample using the WFC3 F160W images fol 
lowing the definition of 

-N-l)\fi\ 



G = 



(3) 



(N-l)Z"\fi\ 

where is the number of pixels within the images and are the fluxes for each pixel sorted in 
ascending order with |/i | < I/2I < ••■ < |//v|, and 



N 



While J2fj<^-^J2\fj\ 



N 



N 



Mr, 



(4) 



(5) 



( Lotz et al. 2004; Abraham et al.ll2007|) ,where Xc,yc is the galaxy' s center and f, are the fluxes for 
each pixel sorted in descending order with |/i | > j f2j > ■■ ■ > | fyl (ILotz et al.ll2004l : [Abraham et al. 



2003). Details of our procedure can be found in A braham et al.l (|2007h . We shows the distribution 
of Gini and M20 for our sample in Figure |6l with randomly selected measurements marked by their 
stamp images. 

Galaxies of different visual morphological types occupy different regions in the Gini-M2o 
diagram (Figures |6]and|7]). Spheroids have higher Gini and lower A/2o, while most irr/mergers have 
lower Gini and high M20. The disks span a wide range in both Gini and M20: those disks with a 
prominent bulge tend to have higher Gini and lower M20 whereas disks with small or no bulges 
have lower Gini and higher M20. 
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For the non-X-ray lEROs, quiescents and dSFGs are well separated in the Gini-Mao diagram: 
nearly all quiescents have Gini > 0.58, while dSFGs have Gini < 0.58 with only a few outliers. 
To further verify this, we re-plot Figure |2]but now with each object color-coded according to the 
Gini coefficient, as shown in Figure [8] This diagram shows that the IR color criteria used in Eq|2] 
to select quiescents also succeeded in choosing objects with Gini > 0.58. The three sources with 
[3.6] - [24] < 0.3 and red K - [3.6] > 1.2 have very low Gini coefficient, suggesting that they are 
not quiescents. 

X-ray lEROs are plotted in Figure[8lb). They show inconsistent classification between their IR 
colors and quantitive morphological parameters. While most of the X-ray sources in the quiescent 
region of Figure |2] show higher Gini ^, nine X-ray sources have Gini > 0.58 yet have [3.6] - [24] 
> 0.3, placing them outside the quiescent region. There are two possible reasons for these sources 
to cross the line: they are quiescents with strong bulges and their 24 fim emission is from central 
AGNs, or they are dSFGs with powerful AGNs appearing in their //-band images and causing an 
increase in Gini. We do notice that, however, most of these sources show a significant bulge except 
one source with ID 1325 appearing to be a point source (Figure [5]). In fact, among all the X-ray 
sources in our sample, 66% (25 sou rces) of them ha ve Lx < 5 x erg s"^ while only 15% (6 



sources) have Lx > 1 x 10 erg s (|Xue et al.ll2011l) . Thus we argue that contamination from the 



AGN in optical bands is not significant for most X-ray sources. Those with high gini are likely 
caused by a significant bulge component. We will further verify this with their rest-frame colors in 
Section [6l 

Though the adopted technique has been used and te sted by many previous works (see, e.g., 
Abraham et al.ll2007 : lcapak et al.ll2007 : Kong et al. 2009), we evaluate the accuracy and reliability 
of Gini and M20 in the following ways. We have already shown that there is good correspondence 
between visual morphological types and their distribution in the Gini - M20 diagram (Figures |6] 
and IT]). Thus the Gini-M2Q system provides a reliable way to automatically distinguish different 
morphological types of galaxies in our sample. Many previous works show that below a certain 
S/N level < S/N >~ 2 , Gini and M^n depend on th e S/N ratio and is useless for distinguishing 
different populations (|Lotz et al.l 120041: lLiskeiil2008|) . Benefiting from the bright cut, the median 
magnitude in the F160W band is // ~ 22.5. And all their images have a mean S/N per pixel < 
S/N >> 2 with a median of<5/A'^>~10 (Lotz et al. 2012, in preparation). Therefore, the Gini 
and M20 for our sample do not suffer S/N effect. As a n independent check, we also calculated 
galaxy central concentration, the concentration index (C. [Abraham et al.lll994h . using also F160W 
images and compared it with Gini (Figure [9]). Our calculated C is closely correlated with Gini, 
consistent with that shown in previous works ([Abraham et al.ll2003|) . This confirms that galaxies 



*The only one that shows very low Gini coefficient, ID 1404, is a high-z source, as discussed in Figure|4] 
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with high Gini indeed are more concentrated. The Gini-C relation shows larger scatter for objects 
with low concentration, mostly dSFGs. This is also expected: measurements of C is based on 
simple aperture photometry, under the assumption that galaxies are circularly symmetric and have 
a well-defined center; however, neither of these assumptions can be fulfilled for the dSFGs in this 
sample. Thus we conclude that our calculated Gini and Mjq provides reliable measurement of 
galaxy concentration and dumpiness for galaxies in our sample. 



5. Photometric Redshifts and Stellar Masses 



5.1. Photometric Redshifts 



Only ~ 23%, 31 galaxies, of the lEROs in our sample have spectroscopic redshifts {zspec), and 
therefore we derive photometric redshifts {Zph) for the rest using the standard SED fitting technique. 
After separating our lEROs into quiescents and dSFGs, we are able to select SEDs accordingly to 
do SED fitting. 



We used the single stellar population (SSP) models of iBruzual & CharlotI (|2003L hereafter 
BC03) to cons truct a template set for the quiescents. We assumed solar metallicity and Calzetti 
extinction law dCalzetti et al.ll2000l) in the range of < Ay < 0.5 mag. For the 10 quiescents with 
spectroscopic redshifts, we obtain very good agreement between Zph and Zspec, with the normalized 
median absolute deviation o"a?mal@= 0.03. Figure [10] shows that the SEDs of the quiescent lEROs 
in our sample are well fit by these SSP models, with best-fit ages ranging from 0.4 Gyr to 4 Gyr. 
Many of them also show a rest-fr ame UV excess, which is likely due to an additional young stellar 
population component (see also lYan et al.l 120041) . However, the UV excess does not affect the 
resulting photometric redshifts, which are mainly determined by the rest-frame optical-to-MIR 
SEDs. 

The dSFGs are much more diverse in that they span a wide range in their [3.6]- [24] colors and 
show a large variety of morphological c haracteristics. We fit these dusty galaxies using the "Easy 
and Accurate Zphot from Yale" (EAzY; [B rammer et al.ll2008l) . We allow for a linear combination 
of six templates which can recover the entire rest-frame color space occupied by observed galaxies. 
An accuracy of (Jnmad = 0-07 was obtained for the 19 dSFGs with spectroscopic redshifts, and their 
SEDs are shown in Figure [TT] 



'Defined as 1.48 x median(\zspec-Zph\/i^+Zspec)) (iBrammer et al.ll2008h 



-13- 



5.2. Stellar Masses 



To obtain stellar masses, we fit the observed B-to-4.5 SEDs o f both the quiescen t and 
dSFGs using "Fitting and Assessment of Synthetic Templates" (FAST: IKriek et alJboOQbh . We 
fixed their redshifts at the derived redshift, and fit the SEDs to a set of simple stellar popu- 
lation models. These models have exponentially declining star-formation histor ies (SFH') with 
0. 1 < r < 10 Gyr and < Ay < 4 assuming Calzetti law. We adopted a Salpeter (ISalpeteilll955b 
initial mass function and solar metallic ity. The stella r masses derived from the best- fitting tem- 
plate are the n scaled to a Kroupa IM F ([Kroupal 120011 ) by dividing the S alpeter stellar ma ss by a 
factor of 1.6 (|Marchesini et al.ll2009|) . We have verified that by using the iMarastonI (|2005|) model 
instead of BC03 models, the stellar mass is systematically smaller with an o ffset of -0.20 and a 
scatter of 0.17. This i s consistent with previous studies for z ~ 2 samples (|Muzzin et al.l l2009l : 



Marchesini et al. 



2009h. We also com pared our stellar mass estimates with that derived using the 



approach in Iwu vts etal.l(l2007[l2008l) . which fitted three different star formation histories: a single 
stellar population (SSP) without dust, a constant star formation (CSF) history with dust (A,, varying 
from to 4), and an exponentially declining star formation history with an e-folding timescale of 
300 Myr and the same range of A,,. The two methods give fairly consistent results with a offset of 
-0.02 and a scatter of 0. 13. 



The mass and redshift distributions are shown in Figure [12] We also derive the mass com- 
pleteness limit as a function of redshift for a 3.6 fim selected ([3.6] < 21.5) sample. This com- 
pl eteness limit is ca l culate d using the full and deeper FIREWORKS catalog following the method 



in 



Marchesini et al.l (|2009l) . While these lEROs span a wide range of 1 < Zph < 4, most of them are 
inl.5<z<2.5 within a narrow peak centered at Zph= 1-8. Their derived stellar masses are more 
striking, with a median value of ~ 1 x IO^'Mq. All this is in broad agreement with the prototype 



lERO sample of lYan et al. These lEROs are clearly among the most massive galaxies. 

At M > 10' (where the our sample is nearly complete; see below), the fraction of quiescents 
increases significantly from ~ 22% at 2 < z < 2.5 to ~ 43% at 1.5 < z < 2 (Table 1), suggesting 
that we are seeing th e quenching process in action in the most massive systems (also see, e.g., 
Cameron etallboil ). 



It is important to understand the completeness of this sar nple in terms of st ellar mass. For 
this purpose, we compare our work to that of FIREWORKS (|Wuyts et al.l l2008b . where the K- 
band selected sample is complete at M > 10'°^A/q to z ~ 3. The mass completeness of a flux- 
limited sample is clearly a function of redshift, and here we concentrate on the redshift range of 
1.5 < z < 2.5 because this is where we are most interested. As stated above, the approach for 
deriving stellar mass for the FIREWORKS galaxies gives consistent results as our estimates. We 
find that our lERO sample has recovered ~ 80% of the FIREWORKS galaxies at M > lO^M© 
within 1.5<z<2.5, as shown in Figure [T3] 
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The other 20% (nine galaxies) have Z85o-[3.6] < 3.25 mag and thus are missed by our se- 
lection. Among these nine galaxies, four of them are X-ray sources (Luo et al. 2008; Xue et al. 
2011), and they also have the highest MIPS 24 fxm fluxes. The SEDs of these four objects have 
strong AGN components, which may have caused their stellar masses to be overestimated. The 
other five galaxies have just slightly bluer zgso - [3.6] colors than the lERO criterion, and only one 
is not detected in MIPS 24 fim. They are clearly at the border of being selected as lEROs. We 
hence argue that the lERO selection using Z85o-[3.6] > 3.25 mag and [3.6] < 21.5 mag is nearly 
complete at > 10^ ^Mq atl.5<z<2.5 except for the bluest galaxies. 



6. Rest-frame Properties of lEROs 



The rest-f r ame c olors permit additional characterization of galaxies. IWuyts et al.l (|2007[) and 



Williams et al.l (12009 ) proposed that quiescents and dSFGs can be separated in rest-frame U -V 



versus V -J: while both quiescents and dSFGs show red U -V colors either due to old stellar 
populations or dust extinction, quiescents are bluer in y-7 due to lack of dust r eddening. We 



comp uted these two colors by interpolating observed data following the algorithm of iRudnick et al. 



(|2003|) . as shown in Figure [141 The rest-frame colors uncertainties were calculated using the input 
photometry through Monte Carlo simulations. For each galaxy, we created mock SEDs using a 
Gaussian distribution and calculated their rest- frame UVJ colors. This procedure was repeated 
1000 times to derive the uncertainties for their color estimation. The average la error ellipse is 
shown in the lower left of Figure [T4b . 

Quiescents and dSFGs in the sample can be well separated in the UVJ diagram (Figure [T4b). 
consistent with our classification based on the observed IR colors and morphological parameters. 
However, there are a few quiescents lying just across the dividing line which may be due to either 
photometric errors or that they suffering extinctions. While we can locate these galaxies in the 
quiescent region by pushing the dividing line to redder V -J colors, many dSFGs will also enter 
the quiescent region. In fact, we verified that the observed K- [3.6] colors for these outliers are also 
redder than other quiescent galaxies at similar redshifts. Interestingly, most of these outliers show 
edge-on disk morphologies in their F160w images, e.g., ID 1 127, ID 2614, and ID 6292 (FigureHj), 
indicating that they suffer more severe dust extinction. X-ray sources can also be separated in 
the UVJ diagram (Figure [T4b): those X-ray sources with high Gini are in the quiescent regions, 
despite that several of them are classified as dSFGs based on the IR colors. We argue that these 
X-ray sources with high Gini may have significant bulges, yet their 24 /xm mainly comes from 
AGN. 



Quiescents and dSFGs can not be separated in the color-mass diagram (CMD) as shown in 
Figure[l5]and Figure[l6] We confine our analysis on their distribution in the CMD to the non-X-ray 
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lEROs in 1.5 < z < 2.5, where our sample is nearly complete for the most massive galaxies. The 
d ashed red line is the red sequence criterion at z = 2 extrapolated from the same criterion at z < 1 



in 



Borch et al.l(|2006|) . Both quiescents and dSFGs in this sample are massive and red. Particularly, 



Figure [T5] shows that the U -V colors span a wide range for quiescents. It is also interesting to 
note that the distribution of spheroids and disks appears similar in the CMD. These quiescent disks 
contribute ~ 30% of the whole quiescent population in 1.5 < z < 2.5^ confirming pervious findings 
dvanderWel etal.boui Icameron etal.lboui iKartaltepe etaklboi ik For the dSFGs, we find that 
disk galaxies tend to be more massive and have redder U -V than the hr/Mergers (Figure [16). In 
fact, as shown in Table [H 70% of the massive dSFGs are disks. These disks are very extended and 
do not show clear signatures of mergers/interactions from visual inspection. M20 for most of these 
disks is smaller than that for Irr/Mergers (Figure 0, suggesting that they are relatively undisturbed. 
These properties are incompatible with the expected compact or highly perturbed morphologies of 
ongoing mergers. Instead, the stochastic accretion of gas rnay play a greater role in t he formation 



of massive disks (IForster Schreiber et all 20061 Dekel et al.ll2009l : Martig et al.ll2009l: lOenzel et al. 



2006L I2OI It iForster Schreiber et al.ll201 1 ). 



7. Conclusions 

With the new Chandra 4 Ms imaging and HST/WPC3 F160W imaging from CANDELS, we 
studied relationships between SEDs, morphologies, and star-formation status for a sample of 133 
IRAC-selected Extremely Red Objects (lEROs) with zssq- [3.6] > 3.25 mag and [3.6] < 21.5 mag 
in the GOODS -S. These criteria select high-mass galaxies at z > 1.5, and our sample is nearly 
mass-complete at > IO^^Mq within 1.5 < z < 2.5. The lEROs show a clear bimodality in their 
IR colors, and we find that they can be classified into quiescent galaxies and dusty star-forming 
galaxies according to their star-formation activities. Such a classification produces very consistent 
results with that using Gini and rest-frame colors (U -V vs. V-J), which indicates that these 
quantities are closely related. Roughly 67% of lEROs in this sample are dusty star-forming galax- 
ies while the remaining 33% are quiescent galaxies. At M > 10' 'M©, the fraction of quiescents 
increases significantly from ~ 22% at 2 < z < 2.5 to ~ 43% at 1.5 < z < 2. Among the whole 
lERO sample ~ 30%, or 38 sources, are AGNs, and the identification of these AGNs are shown 
to be very important in understanding the discrepancies among different classification methods. A 
number of galaxies are identified as star- forming galaxies based on their 24 fxm emission, yet are 
in the quiescent region of the UVJ diagram; we show that most of them are actually AGNs. 

The morphological study in this work provides important clues to massive galaxy evolution 
at z ~ 2. The quiescent and dSFGs in this sample show distinctive morphologies. While all quies- 
cents have a prominent bulge, most of the dSFGs are either disks or irregular/mergers with no or 
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small bulges. This suggests that a prominent bulge is necessary to form a massive quiescent galaxy. 
We also show that the two populations can be well separated based on the Gini coefficient. The 
quiescents have higher Gini coefficients than the dSFGs, indicating that the formers are more con- 
centrated. We argue that the quenching process for star formation must lead to or is accompanied 
by the increase of galaxy concentration. 

Roughly 30% of the quiescents also show an extended disk component, confirming recent 
findings with a larger sample. Moreover, we also find that ~ 70% of massive dSFGs with > 
IO^^Mq at 1.5 < z < 2.5 are classified as disk galaxies. Most of these disks appear relatively 
undisturbed with no clear signatures of mergers/interactions, based on both visual inspection and 
quantitative measurements. The prevalence of disks among massive galaxies at z ~ 2 challenges the 
merging scenario as the main mode of massive galaxy formation, and suggests that the stochastic 
accretion of gas plays a greater role. 

This work is based on observations taken by the CANDELS Multi-Cycle Treasury Program 
with the NASA/ESA HST, which is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS5-26555. We gratefully acknowledge R. Abraham for 
access to his morphology analysis code. This work is also supported under the National Natural 
Science Foundation of China under grants (10878010, 10221001, 10873012 and 10633040) and 
the National Basic Research Program (973 programe no. 2007CB815404 and 2007CB815405). 
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Fig. 1. — Observed-frame zg5o - [3.6] vs. [3.6] - [4.5] color-color diagram for all IRAC- selected 
([3.6] < 21.5) ga laxies in the GQQDS -S. Observed-frame color tracks for nonev olving E/SO, Scd 
SEDs (CWW E icolananeLalJ[l98Ql). and a dusty star-forming galaxy rHRlO. IPey et al.lll999l : 



Stem et al.ll2006h SED at < z < 4 are plotted using color-coded lines. The lERO color definition, 



Z85o~[3-6] > 3.25, is marked with a black dashed line. The red Z850~ [3-6] color can be produced 
by both old and dusty galaxies at z > 1 . 



Table 1. Visual classification results for lEROs in our sample. 



Type 


All 


Spheroids 


Disks 


Irr/Mergers 


Total 


133 (38)^ 


34 (8) 


59(17) 


40(13) 


Quiescents 


43 (8) 


25 (2) 


15(4) 


3(2) 


dSFGs 


90 (30) 


9(6) 


44(13) 


37(11) 


Massive'' Quiescents at 1.5 < z < 2.0 


10(6) 


4(2) 


5(3) 


1(1) 


Massive Quiescents at 2.0 < z < 2.5 


4(0) 


2 


1 


1 


Massive dSFGs at 1.5 < z < 2.0 


13(5) 


1(1) 


8(3) 


4(1) 


Massive dSFGs at 2.0 < z < 2.5 


14 





11 


3(1) 



^The number of X-ray sources. 

"'M* > IO'^A/q, where our sample is nearly complete. 
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Fig. 2. — Tlie [3.6] - [24] vs ^-[3.6] color-color diagram for galaxies in our sample, which is 
divided into non-X-ray (left panel) and X-ray sources (right panel). Upper limits correspond to 
galaxies fainter than 1 1 /iJy {3a limiting flux density) in the 24 /xm image. The color tracks for 
the CWW E/SO and HR 10 SEDs at 1.5 < z < 3 are shown in red and blue lines. The blue and red 
filled circle indicate the locus of their colors at z = 1 .5 and z = 3, respectively. The [3.6] - [24] color 
provides an estimate of SSFR while the ^-[3.6] color is sensitive to dust extinction. The [3.6] 
- [24] color distribution is clearly bimodal. All the 24 //m- undetected lEROs are at [3.6] - [24] 
< 0.3, most of which also have bluer ^- [3.6] colors. We propose a preliminary color criteria of 
[3.6] - [24] = 0.3 and ^- [3.6] = 1.2 to separate quiescent and star-forming galaxies. The central 
AGN can also have significant contribution to their 24 fim emission, and thus this criterion may 
not apply to some X-ray lEROs. 



Spheroids Disks Irregulars/IVIergers 




Fig. 3. — For the visual classifications, we simply divide the sample into three broad morphologi- 
cal types: spheroid on the left, disk morphologies in the middle, and irregular/mergers on the right 
panel. 
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Fig. 4. — HSTfWFCS F160W images for the quiescent galaxies (quiescents) in our sample. Qui- 
escents are defined by Eq 1 and lie within the dashed square in Figure |2] The source IDs (from 
FIREWORKS catalog), M20, and Gini are labeled in each panel. X-ray sources are marked with an 
"X" in the bottom-left comer. Both non-X-ray and X-ray sources are sorted in order of increasing 
^-[3.6] colors. The size of each postage map is 4" x 4", and 1" corresponds to ~ 8.5 kpc at 
z ~ 2. Nearly all objects have a prominent bulge component with ~ 30% of them also showing 
an extended disk. The only outlier is Object 1404, which shows a very diffuse morphology with 
a photometric redshift z ~ 3. At this redshift, the MIPS 24 fim probes the rest-frame 6 jim, and 
even star-forming galaxies has low brightness. We thus argue Object 1404 is not a z ~ 2 quiescent 
galaxy. 
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Fig. 5. — HSTfW¥C3 F160W images for the dusty star-forming galaxies (dSFGs) in our sample. 
The dSFGs are all galaxies that lie outside the dashed square in Figured X-ray sources are marked 
with green "X". Fireworks ID, M20, and Gini are shown in each postage stamp. Both non-X-ray 
and X-ray sources are sorted in order of increasing K-[3.6] colors. Most of these dSFGs have 
disk or irregular/merger morphologies, with only a few of them, mainly X-ray sources, showing a 
strong bulge component. 
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Fig. 6. — Gini vs M20 for the galaxies in our sample. Quiescents and dSFGs are defined based on 
Eq [2l Non-X-ray sources are denoted with circles and triangles while X-ray sources are denoted 
with filled circles and triangles, respectively. The quiescent and dSFGs are well separated in this 
diagram (upper left insert). The postage-stamp images for a randomly selected subsample represent 
morphologies for galaxies in various location of this diagram. The quiescents have high Gini and 
low M20, and the dSFGs primarily have low Gini but cover a wide range of M20. The quiescent and 
dusty star forming galaxies with no X-ray detection are well separated by Gini = 0.58 with only a 
few outliers. The X-ray sources do not show such consistent classifications between morphologies 
and their IR colors. 
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Fig. 7. — The distribution of all the lEROs with different visual morphological types in the Gini- 
M20 diagram. Quiescents are denoted with red colors while dSFGs are in blue colors. Their 
distribution shows great agreement with the visual classification results. Spheroids have higher 
Gini and lower M20 while Irr/Mergers have lower Gini and higher A/2o. Interestingly, quiescent 
disks show primarily higher Gini than the star-forming disks, consistent with that the former have 
a prominent bulge. 
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Fig. 8. — [3.6] - [24] vs .^-[3.6] color-color diagram for the lERO sample with galaxies color- 
coded according to the Gini coefficient. A nearly perfect correlation exists between [3.6] - [24] 
color and Gini. Such a correlation, does not exist for the X-ray sources in the right panel. 
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Fig. 9. — Gini coefficient (Gini) vs. central concentration (C) for lEROs in our sample (points as 
in Figure |7J. The solid line corresponds to unity slope, while the dashed lines are ± 10% offset 
relative to the solid line. 



-30- 



100 

1 

0.01 



age=0.5 



id=317 

zp=1.42 

zs=1.34 



age=0.4 



id=1127 
zp=1.40 



age=0.4 



id=2614 
zp=1.56 



age=3.0 



id=4207 
zp=1.48 
zs=1.98 



age=0.5 



id=4982 
zp=1.66 
zs=1.61 




100 








age=0.5 




1 


■ ^ -jf^ 


id=6097 ■ 






zp=2.04 


0.01 




zs=1.90 




0.5 1.0 


5.0 



age=0.3 






id=6187 




zp=1.58 




zs=1.61 


0.5 1.0 


5.0 



0.5 1.0 



5.0 



Fig. 10. — SED fitting to measure photometric redshifts for the quiescents with SSP models as 
templates. The source ID, age of the best-fit SSP model , photometric redshift (zph) and spectro- 
scopic redshifts(Zipj,c, if available) are labeled in each panel. Most of these quiescents can be well 
fit with SSP model templates with ages ranging from 0.4 Gyr to the age of the Universe of 4 Gyr 
at that redshift. "X"s denote the X-ray sources in Figure 5] 
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Fig. 11. — Example SEDs for the dSFGs with EAzY spectral fits overlaid. Here we only show the 
SEDs for those dSFGs with spectroscopic redshifts. Data points with la errors are shown in black. 
The red line shows the linear combination of template SEDs that best fit the data. The source ID, 
Zph and Zspec are labeled in each panel. 
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Fig. 12. — Redshift and stellar mass histograms of the full sample. The dashed line corresponds to 
the mass completeness limit for a 3.6 /zm selected sample ([3.6] < 21.5) at 1 < z < 4. The redshift 
range of the lERO sample spans from 1 to 4 with most of them located at 1.5 < z < 2.5, indicating 
that the simple color criterion of z^so - [3.6] > 3.25 can effectively select massive galaxies at z ~ 2. 
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Fig. 13. — Z85()- [3.6] vs stellar mass for all galaxies at 1 .5 < z < 2.5 in the FIREWORKS catalog. 
Red points are MIPS sources. The purple open squares are lEROs in our sample (z85o- [3-6] > 3.25 
and [3.6] < 21.5). There are 49 galaxies with > 1 x lO'^M©, which all have [3.6] < 21.5. 
About 80% of them (39) are selected as lEROs with Zs5n ~ [3.6] > 3.25. Among the nine galaxies 
not selected as lEROs, four of them are X-ray sources (|Luo et al.ll2008l : IXue et al.ll201 lb . which 
also have the highest MIPS 24 fim flux and bluest Z850 ~[3.6] colors: Z850 ~[3.6] < 3.1. Their 
SEDs have strong AGN components which may cause their stellar mass to be overestimated. The 
remaining five galaxies just missed being lEROs according to Z850~ [3.6] being too blue. Only one 
of them is not detected at 24 /im. 
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Fig. 14. — Rest-frame U — V vs. y-7 color-color diagram for non-X-ray (left) and X-ray (right) 
lEROs in our sample. The error ellipse is shown in the lower left of the left papel. Objects are 
color-coded based on Gini coefficient. Quiescent and dSFGs can be well separated in this diagram 
with the dashed line. The locations of high- vs. low-Gini objects on opposite side to the lines 
confirms that morphology and rest-frame color-color yield generally consistent classifications for 
both populations. The distributions overlap somewhat, but mostly near the boundary. Both massive 
quiescent and dusty galaxies have red U -Y colors, but the quiescents are bluer in y-7, implying 
they are not suffering significant dust extinction. X-ray sources with different Gini values are about 
as well separated as non-X-ray sources (right panel). 
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Fig. 15. — Rest-frame U-V color vs stellar mass for quiescents in 1 .5 < z < 2.5 in the sample. HSTfWFC3 
F160W images are also shown. Visual classification results for these quiescents are shown in the upper left 
p anel. The dashed re d line is the red sequence criterion at z = 2 extrapolated from the same criterion at z < 1 
in borchet all ( 2006 ). 
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Fig. 16. — Rest-frame U-V color vs stellar mass for dSFGs at in 1.5 < z < 2.5 in the sample, also 
shown are HSTfWFC3 F160W images for a randomly selected dSFGs sample. Visual classification results 
for the dSFGs are shown in the upper left panel. Most of the dSFGs with > 1 x lO^^M© show disk 
morphologies. 



